E A 


DE E o co Si 


rias 


-.. 


Es 


a” “mM 


VARÕES DE AÇO DE ALTA RESISTÊNCIA 


PARA BETÃO ARMADO TIPO A. 40 
(«Heliaço» e «SNT 40») 


TENSÕES DE SEGURANÇA 
Ra = 2400 kg/cm? 
e 
Ra = 4000 kg/cm? (SISMOS) 


Literatura e Tabelas de Cálculo 
fornecidas gratuitamente 


HELIAÇO, LDA. 


R. DE S. JULIÃO, 41, 4.º D, — LISBOA 


Tel.: 36 7817-3237 06-369477 


DAIÃO ELÉCTRICA PORTUGUESA 


S. A Ri bl. 


SEDE 
Rua de Alexandre Herculano, 244 — PORTO 
TELEFONES: 27676/77/78/19— Est. 90 


DELEGAÇÃO 
Rua Rosa Araújo, 35 — LISBOA 
TEL.: 7311 81/82/83/84/85/86 - 733333 - 48491 


ELECTRICIDADE DO LINDOSO 
E DAS CENTRAIS DO FREIXO 
E DA CACHOFARRA 


“ 
A UNIÃO ELÉCTRICA PORTUGUESA distribui, 


pela mais extensa rede de alta tensão em Por- 
tugal (cerca de 2500 km), energia eléctrica em 
70 concelhos, nos distritos de VIANA DO CAS- 
TELO, BRAGA, PORTO, AVEIRO, COIMBRA, 
VISEU, LEIRIA, SETÚBAL e ÉVORA. 


º “e... . +. 

nene." * “0a. “eau... 
PO EMO DO te e e o e e o e o o 0 0 000 0 66 
COCA e A a 0 4 ... 
COCO e q 0 é 


MiKMOEXAO TX Md mM ri SZZEANErAKHXAMAX ANA, 
AXEL TA AX A O A KANAAN ZLAHXAHA. 
MAXX x 7 x A IMZZHAXZXAKAKAMAAAARA, 
27000: 00600-0DD Cc 04009: 90,90:0:0 9:0 0,0:0,0:0 
ARA ZXAX A x x 0 AESA 
92 00,000070 0 O “COD 0:;0,0-:0/0-0:6 3:00:00 
DX É a A SR EEE IXIEAMARARDA. 
ATA TV ESC a 1 PE AD ND PRE A e E Sr TC AMA AX. 
CCC Ci. “RCA MMMAÕÃ ÃO..." + 
ecc. E Cad a EM AOÕJÃÃjÃ££te........ 
cce ARITNANANAN TT TaUadHT A OõCÕJÃãÃãÃãÃÃrere. en..." 
coco E FF Pd Ni A DUDE oe... ... 
AH Wii cats Vaca da A Dx XANAX 
aAaMTXEM ESSES A FI E E CZK AM AA, 


+ 
= 00 
> o, .*. 


MAQUINAS E :: 
EQUIPAMENTOS  :: 
AMERICANOS ::: 


. . - 
eco. Was 
000000 War 
“ee... 


ui. Telelelelels MEPFREDENTANTES ENCLUuBvoOs:. 
cet AUTO-LUSITANIA (arreoo DUARTE, LOM.) 
etetete teto teto et eta t ato Avenida de Liberdade, 73-1º e 2.º, 77 — LISBOA 2 


e..." 
coco ce ce e e o a 0%. 


e... .......“ 


0200000 0900000006 Tema 00 00 0:0:0:0:0:0:0:0 
esccce ce ce ce e e o ee. O at CC o 0 0 0 0 0 00.0 
ecc. 0 «UC CCC CCC 0 0 0. 4 
7000000000 000 000000 0000.» 90000 00000000 0:00 0 000.0 
“ecc... 00 «MC 0 0... | 
CCC. e nn 0 q 
escore co o o ve ce e co e a e oo ce e e es o e o un 0 0 0 0 0 0 0 0 0 0 0 0,00, 0,0,0.0,6 
0 0 0 0 0 0 e 0 e 00 00 000000000 60 0 00 0 0 0 0 0 0 0 0 0 6 00 0 0 00 0 0 o 0 00” 
cocos o co cv cc e o e ce ce e e e o e e o e vo o ne e a ce o o 0 0 O 0 0 0 0 0 0 0 0.4 


COCA A e A O A Un 0 0 0 0 6 4 
orcs co sc e e e e e ne e nn ue e nn cu uu ou ou ue o o O o DO 0 0 0 0 0 4 6 


LIMENTO SELL 


ES 


E 


Fábrica: 
Outão —- Setúbal 
Telef. 228824/5 
Escritório: 
KR, do Comércio, 56-3.º 
LISBOA 2 
Telef. 328201/2/3 


NUMERO 3d JANEIRO DE 1968 


Director 
Vitor M. M. Pinto 


Administrador 


Vitor Coias e Silva 


Redactores 


Ana Maria Sobral 
A. J. Tavares Pires 
J. C. Alves da Silva 


João M. Pereira 


Editor e Proprietário 


Associação dos Estudantes 


do Instituto Superior Técnico 


Redacção e Administração 


Avenida Rovisco Pais 


Telefone 77 5449 


Lisboa 1 
Assinaturas 

5 N.os 10 N.os 
Continente e ilhas — 70500 150300 
Ultramar, Brasil 
e Espanha — 80500 140$00 
Estrangeiro — — M$00 1505800 


Número avulso — 20800 


Não se publica em Agosto 
e Setembro. 

Os artigos assinados são da 
exclusiva responsabilidade dos 


seus autores. 


Composto e impresso na Jipografiu 
forge Fernandes, Lda. — Hua da 
Cruz dos Poiais 103 — Lisboa 


TECNICA N.º s74 


SUMÁRIO 


189 


201 


213 


219 


231 


241 


239 


247 


CAPA; 


ANO XLII VOLUME XAA PUBLICAÇÃO MENSAL 


CARLOS PORTELA — Análise de alguns efeitos de 
perturbações de tensão nas máquinas sincronas 


Analysis of the effects of voltage disturbances over synchro- 
nous machines 


ARTUR RAVARA — Aplicação dos computadores ao 
cálculo dinâmico de estruturas — Application of com- 


puters in the dynamic analysis of structures 


ALEXANDRE CERVEIRA — Distância mínima nos có- 
digos de Bose-Ghaudhuri-Hocquenghem. Método 
para a determinação de um limiar 


Minimum distance for Bose-Chaudhuri-Hocquenghem's codes 
A method for the determination of a threshold 


LEOPOLDO CORTEZ — Princípios gerais e métodos 
de investigação na flutuação de minérios 


Ore flotation: general principles and investigation methods 


PAUL DUPONT — Procédé suprimant |'emploi clas- 
sique des cintres dans le construction des arcs 
et voutes de ponts en béton armé ou non armé 
Procedure allowing the suppression of centering in the cons- 
truction of concrete or reinforced concrete bridge arches and 
vaults 


Do Mundo Técnico 


A adaptação da formação dos engenheiros às exigências do 
crescimento económico — Nota sobre a reparação de fendas 
em chapas de grande espessura 


Notas Informativas 


Elementos sobre a produção e o consumo de energia na rede 
eléctrica nacional 


al - 
Bibliografia 
Torre de carregamento de navios. Débito 750 t/hora,. Em funcionamento no 


terminal de açúcar em Lourenço Marques — Frojectada e construida em 


Alverca por Construções Metalomecânicas MAGUE, S. A, R.L, 


187 


ÍNDICE ALFABÉTICO DOS ANUNCIANTES 


Pág, 
REG suma a ga og a XVII 
A. Johnson & Co (Portugal), Ltd. : XXH 
Aguiar e Melo, Lda. . . +... o XXI 
Asea Eléctrica . . V, VII, XXV e XXVII 
Auto-Lusitânia - 2.º da capa 
Cerâmica de Vale de Lobos, Lda. . XVI 
Cimento Secil. . . .. cc... « 2.º da capa 
Companhia de Cimentos Brancos . VI 
Construções Técnicas . . . . .... XHI 
Dida smos sd Ea, SE Gê E XVI 
Electrotécnicos Reunidos . « X 
Electro-ÁArco, Ltdd . «wo. +. “ H 
Empresa de Cimentos de Leiria . . XXVI 
Empresa Electro Cerâmica . . +. . |. XXV 
Empresa de Produtos Asfálticos . . VIII 
Empresa de Sondagens e Fundações 
Teixeira Duarte, Lda. . . . . . .3dacapa 
Empresa Técnica de Equipamentos 
ENO é Es sé gg gs sa IX 
Empresa Ultramarina de otida gas 
e EUNdAÇÕES E mé cpm s à XIV 
Fanafel . . +... E RUE ME E 3 XIV 
Fundações Franki, Hã. E EMA Á I 
Guedes & Almeida, Lda. . . ... XVI 


188 


Pág 
Heliaço 2.º da capa 
LE meto Nm nas a XVI 
Jayme da Costa, Lda. . . . .... XV 
Ita Rule; ss GEES a V 
J. F. de Azevedo e Silva +... IV 
; E Sivadoe lmdo o ss vigas q XXVI 
MAGUE sta sc E A XI 
Novobrá, Lda. « cus sumo q XVI 
Ma o as ums & Em a XH 
SIMÓN ss sema * és - 4.º da capa 
Soc. Industrial Metalúrgica . . .. VII 
Sociedade Portuguesa Cavan . . +... VI 
Sondagens e Fundações A. Cavaco . H 
Sondagens Ródio, Lda. . . +... VII 
DUGRENHE ms ese nm ams » H 
União Eléctrica Portuguesa - 2.º da capa 
Wild Portugal, Lda. +... . +... HI 


TECNICA N: Sid 


qui st Es EM 


ia 


. 
Ea 


= » 


aqu eotigas 


k 
2 


- 


TE E CANETA 
po Paes Si A a E 


Mo e 
= e º 
Ja 


” " a 


Companhia Cimento Tejo — Alhandra 


Assente sobre Estacas FRANKI 


“UNDAÇÕES FRANKIL?* 


RUA BRAAMCAMP, 13-3.º — LISBOA 1 — TEL. P, P. C. 5341 12/3/4 — ENDEREÇO TELEGRÁFICO: FRANKI — LISBOA 
PRÉDIO SANTOS GIL, 6.º ANDAR — SALA 10 — TEL, 91167 —- CAIXA POSTAL 1397 


ENDEREÇO TELEGRÁFICO;, FRANKI-LOURENÇO MARQUES 


TECNICA | 


SOPECATE 


y SONDAGENS 
Ê FUNDAÇÕES 


Escritório—R, DO ARSENAL, 146, 2.º Esq. Dt.º e Fr.'º 
Estaleiro e Oficina — RUA PEREIRA HENRIQUES, 58 


= 


MERDA TaTrEras 


DA 
1 


SERVIÇOS ADMINISTRATIVOS 34010 
TELEFONES + SERVIÇOS TÉCNICOS 320208 


ESTALEIRO E OFICINAS 381238 
LISBOA 


ELECTRO-ARGO, L.” 


ELÉCTRODOS + POSTOS + ACESSÓRIOS 
PARA A SOLDADURA ELÉCTRICA 


e ESTUDOS, PESQUISAS E CAPTAÇÕES DE AGUA 


e SONDAGENS 
= 


MATERIAL APROVADO PELOS 


LLOYD'S REGISTER OF SHIPPING 
BUREAU VERITAS 
GERMANISCHER LLOYD 


* REBAIXAMENTO DE NÍVEIS AQUÍFEROS 


e FUNDAÇÕES 


LISBOA 
Rua Silva Carvalho, 239 
Tels. 6836 49-684893 


FÁBRICA E LABORATÓRIOS 
Venda Nova — Amadora 
Tels. 9706 13-971733 


e) SONDAGENS E FUNDAÇÕES A. CAVACO. LDA. PORTO 
R. RODRIGO DA FONSECA. 62, RC ESQ. — LISBOA -1 Rua do Bolhão, 216 
- TELEF. 5611 71/72/73 Tel. 21277 - 394473 


TECNICA | 


O novo Nível automático 
de engenheiro 


WILD NAK 2 


com circulo horizontal 


Para nivelamentos de média e grande precisão; 
com micrómetro de lâmina de faces plano- 
-paralelas para nivelamentos de grande precisão. 


Erro médio para 1 Km de 
nivelamento duplo: + 1,5 mm 


Erro médio para 1 Km de 
nivelamento duplo com micró- 
metro: + 0,4 mm 


Luneta de 30x de aumento 
com imagem direita de obser- 
vação 


Solicitem catálogo ou demonstração 
aos representantes exclusivos 


WILD PORTUGAL, LDA. ciimaon ini 


TECNICA III 


SAUNIER 
DUVAaAL 


ILUMINAÇÃO LAMPADAS LÂMPADAS 


PÚBLICA INCANDESCENTES ESPECIAIS 


E FLUORESCENTES 


-SIMEL 


TERMINAIS ARMADURAS 
E LIGADORES TIPO OLHO DE BO! 


TAYLOR-" 
HOBSON 


BRAÇADEIRAS 


APARELHAGEM 
DE METROLOGIA 


METROLOGIA E MICRÔMETROS 
PARA A INDUSTRIA 


PROJECTOS 


ESTUDO E CONCEPÇÃO 
DE INSTALAÇÕES ELÉCTRICAS 


ILUMINAÇÃO 
DECORATIVA 
E INDUSTRIAL 


CONTACTORES CANALIZAÇÕES 
DISJUNTORES ELECTRICAS 


FUSIVEIS DE APC 
E INTERRUPTORES 


MATERIAL 
DE GRAVAÇÃO 
EM FITA PLASTICA 


ANTENAS FECHADURAS 


HELLERMAN 


PROJECTORES 
FITAS DE AFILACAR FITAS ADESIVAS ÓPTICOS E 
MICRO - INDICADORES 


TELEFONES 
DE PORTA 


MATERIAL DE SOM 


PRODUÇÃO 
CONCEPÇÃO E FABRICO DE QUADROS, 
CONVECTORES E EQUIPAMENTOS 
PARA INSTALAÇÕES ELÉCTRICAS 


MONTAGENS 


EXECUÇÃO E MONTAGEM 
DE INSTALAÇÕES ELÉCTRICAS 


JE DEAZEVEDO E SIEYA. &S& CG: E! 


TECNICA IV 


Rua S. Francisco Sales, 2 Telef. 654165 - 


Lisboa-1 


pás 550 esc. 


CARACTERÍSTICAS: Tubo de aço, articulado em qualquer posição, várias Incli- 
nações, fbclimente trensportável, não ocupe espeço quando 
fochado, prenchets de 1,20 X 0,80 de medeiro, em tola. 
Pintura de estufo, 


Em exposição, outros modelos de madeiras e metólicos e todo o mete- 
rial de consumo e equipamento de desenho, escritório e papelsris 


Distribuidores exclusivos, em Lisbos, do estirador da Fábrica Olaio 


J. J. SILVA & IRMÃO, LDA. 


rus silvo carvalho 232-b » lisbos - telefons 652750 | 680671 


- 
: 
“a 
(8) 
= 
« 
= 


TECNICA 


O EMPRECO DO CIMENTO BRANCO... 


permite acabamentos 

mais perfeitos, mais 

duradouros e muito 
mais económicos 


Estude a vantagem do emprego do 


CIMENTO BRANCO LUSO 


IBRA 


Consulte os distribuldores gerais 


No Sul: SCIAL —T. do Corpo Santo, 15 — Telef. 66 61 86 — Lisboa 
No Norte: SCIAL — R. do Bonjardim, 205 — Telef. 25779 — Porto 


POSTES DE BETÃO PARA 
LINHAS ELÉCTRICAS 
MOSAICOS E AZULEJOS 


COLUNAS PARA 
ILUMINAÇÃO PÚBLICA 


MANILHAS PARA ESGOTOS 


das! 
as E e, 
E ES r É] 
: 
4 
+ 
ms 
b 
. 4 
E o a! 
i 
a h 


RUA D. ESTEFÂNIA, 94-A-LISBO AA -TELEF. 47812-501 


TECNICA MI 


b 

| 

| 
pes 
' 4 ] 

] 


29 


SONDAGENS RÓDIO, L.* 


LISBOA 


RUA DE S. BENTO, 644-3.º 
Telefones: 68 80 96/7/8 4 Telegramas: SETANSOL 


SONDAGENS GEOLÓGICAS 

ESTUDOS GEOTÉCNICOS 

BS TA CAS GUNATA 

CAPTAÇÃO DE ÁGUAS 

CONGELAÇÃO DO TERRENO 
IMPERMEABILIZAÇÃO DE OBRAS HIDRÁULICAS 

CONSOLIDAÇÃO DE FUNDAÇÕES 
ABAIXAMENTO DO LENÇOL DE ÁGUA 


ie / | — 
Engenheiro Consultor: Walter Weyermann 


Sociedad InoustRiaL Merauncica 


Responsabilidade Limitada 
(REGISTADO) 


SERRALHARIAS, 


CALDEIRARIA, 
FERRARIA, 
FUNDIÇÕES 


um 
ESCRITÓRIO 


Rua de 5. Tiago, 13 


LISBOA 


Telefone 86 19792 


TECNICA VII 


Celulite EE 


CONSTRUÇÕES 


E É PAVIMENTOS 
Nida = FF TECTOS 
E qa 
Cy am E 
ao — — a PAREDES 
== E a 
sm — CÂMARAS 
= dm: 7 FRIGORÍFICAS 
= * 
ISOLAMENTOS 
INDUSTRIAIS 


Terraços e pavimentos isolados com CELULITE no 
edifício da Comp.' de Seguros Fidelidade — P. das Águas Livres — Lisboa 


SETH, LDA. + RUA FILIPE FOLQUE, 10-1.º + TEL. 5301 56 — LISBOA 


| |: 1h DE PAn 
Para impermeabilizar terraços, E estica dos 
paredes, fundações, etc. <q 


IMPERMEABILISADOR 


COLADEM DE 
TACOS DE MADEIRA 


TELEF. 530156  TELEG. EPALDA — LISBOA 


PABRICA EM SACAV EN 
TECNICA VIII 


HEWLETT PACKARD 


aparelhagem electrónica 
de medida e registo 


Osciloscópios, Oscilógrafos, Voltimetros 
electrónicos, Geradores de sinal, 
Contadores electrónicos, Registadores 
digitais, Fontes de alimentação 


Oscilocópio HP 


Contador Electrónico HP 


Comunicações 


equipamentos 


e sistemas completos 
+ Cuidadosamente 
projectados 


+ Amplamente 
experimentados 


* Completamente 
integrados 


THE TECHNICAL MATERIEL CORPORATION 


Registador Sanborn 


Registador digital fura drdi dd 
Si O dei dei 
e. aaa a ao 6 


JELEG NSHRA EMPRESA TECNICA DE EQUIPAMENTOS ELEGTRIGOS, SARL. 


TECNICA IX 


TELEFONBAU UND NORMALZEIT 


TECNICA X 


Os P. P. C. da Telefonbau são fabricados em 3 
modelos com capacidades de 30 a 200 exten- 
sões e permitem a ligação de serviços especiais 
para a indústria hoteleira, nomeadamente serviço 
de despertar automático, chamada de criados, 
bloqueamento da linha de rede, inter-ligação 


com central automática privativa. 


ef RTP re RNNAME 
EM ENG: 
DETIDOS | “a 
UU : IUJN] IUAUR 2, [ da. 


NUMERO 374 JANEIRO DE 1968 


ANO XLIII 
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C.D.0, 624.343.32.015 


ANÁLISE DE ALGUNS EFEITOS DE PERTURBAÇÕES 


DE TENSÃO NAS MÁQUINAS 


RESUMO 


Neste artigo Jag-se a análise de alguns efeitos das va- 
riações bruscas de tensão nas máquinas sincronas, pro- 
curando bisicamente um tratamento simplificado, para apli- 
cação numérica sistemática, 

Analisam-se, designadamente, os aspectos de estabilidade 
dinâmica e de corrente nos enrolamentos induzidos, 

Em condições em que sejam aproximadamente válidas 
as hipóteses indicadas, os resultados numéricos apresenta- 
dos permitem averiguar as condições de estabilidade de má- 
quinas siíncronas perante uma perturbação de tensão e as 
consequências dessa perturbação sob o ponto de vista da 
sobreintensidade no induzido, 


0. INTRODUÇÃO 


Neste artigo faz-se a análise de alguns efeitos 
das variações bruscas de tensão nas máquinas 
síncronas, procurando basicamente um tratamento 
simplificado para aplicação numérica sistemática, 

Analisam-se, designadamente, os aspectos de 
estabilidade dinâmica e de corrente nos enrola- 
mentos induzidos. 


1. HIPÓTESES BÁSICAS 

A análise numérica baseia-se essencialmente 
nas seguintes hipóteses : 

a) Quanto à perturbação de tensão. 


Supõe-se que a tensão anterior, equilibrada, 
tinha um valor V, e que se dá uma pertur- 
bação brusca, continuando a tensão equili- 
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SÍNCRONAS 


por CARLOS DE MEDEIROS PORTELA 


Doutor em Engenharia Electrotécnica 
Chefe do Centro de Cálculo da Companhia 
Nacional de Electricidade 


SYNOPSIS 


This paper analyses some effects of sudden voltage dis- 
turbances over synchronous machines, basicaly trying a sim- 
plifled treatment for the systematic numerical application. 

Analysis is performed, namely, on the aspects of dyna- 
mic stability and of current in armature winding. 

Whenever the hypotheses considered are approximately 
valid, the numeric resulta hereby presented allow to check 
the conditions of stability in synchronous machines, within 
a voltage disturbance, and the consequences of that distur- 
bance, in viem of the over-intensity in lhe armature, 


brada (componentes inversa e homopolar 
nulas) e passando a componente directa a 
um valor V, sem variação de fase nem de 
frequência, mantido durante um tempo T, ao 
fim do qual a tensão retoma o valor inicial. 


b) Quanto ao comportamento da máquina síncrona. 


Supõe-se que o comportamento dinâmico da 
máquina síncrona, trifásica, é tradutível por 
uma inércia constante, uma potência mecá- 
nica constante, uma f. e. m. «equivalente» de 
amplitude constante e fase rigidamente ligada 
à posição mecânica do rotor e uma reactância 
«equivalente» constante. 


A generalidade dos resultados obtidos é apli- 


cável quer ao caso do funcionamento como gera- 
dor quer ao do funcionamento como motor. 
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O tratamento das grandezas associadas à com- 
ponente inversa da tensão durante a perturba- 
ção, que não é aqui analisado, pode, sob diver- 
sos aspectos, ser feitô separadamente, quer 
quanto a estabilidade quer quanto às perdas 
globais. 


2. RELAÇÕES ADOPTADAS PARA O PRO- 
BLEMA DA ESTABILIDADE 


Sejam, dentro das estilizações anteriores: 


t o tempo (em segundo), a contar do 
início da perturbação de tensão 

T | a duração da perturbação de tensão 
(em segundo) 

WD. a pulsação eléctrica, igual a 27f, sendo 
f a frequência em Hz 

6) a velocidade angular do sistema rota- 

tivo (em rad seg”) 

o «ângulo de potência» da máquina 

(em rad), supondo 2 positivo, no caso 

dos geradores para f.e.m. em avanço 

em relação à tensão, e no caso dos 
motores para a f.eem. em atraso em 
relação à tensão 

o valor de 2 antes da perturbação de 

tensão em regime estabilizado 

E a f.e.m. «equivalente» referida à ten- 
são entre fases (em Volt) 

V. a tensão aos terminais fora do período 
de perturbação da tensão, referida à 
tensão entre fases (em Volt) 

V a tensão aos terminais (em Volt) 


4 


e 


|, a corrente no induzido antes da per- 
turbação 
I a corrente no induzido no caso de 


haver perturbação 
X a reactância equivalente da máquina, 


por fase 
W a energia cinética do sistema rotativo 
(em Joule) 
PR 4 
Va 
sê. 
Vê 
Es o VE 
2WX' 
Tem-se que, durante a perturbação da tensão, é 
d? 3 À 
1) = A (sin 9 —v sind) 
dt? : 


“ — 


dj et nica 
2) a TV2A |] (sin d, — vsin 9) dô, 
Jo 
e após terminada a perturbação de tensão é 


3) — = A (sind — sin 9) 


dj [ fT o * . * s 
4) Ei === V24 1) (sin %0—vsin 9) do + 
t 


6) 


= 
+[ (ain ds — sin) dê |. 
x 


3. ANÁLISE DAS CONDIÇÕES DE ESTABI- 
LIDADE DINÂMICA 


Como se verifica fácilmente, a condição da 
tensão V poder ser mantida indefinidamente sem 
instabilidade corresponde à condição de a área a 
tracejado (algébrica) indicada na fig. 1 ser posi- 
tiva. 


É (grau) 


Y 
á, SH rod) 


Fig. 1 — Relação entre o ângulo de potência e a potência 

trocada com a rede, P, antes da perturbação de tensão 

(v==1), e durante a perturbação de tensão (v = v). Repre- 

sentação correspondente à verificação de estabilidade para 
uma perturbação de duração infinita. 


O valor mínimo da tensão V que pode ser 
mantida indefinidamente após a perturbação sem 
perda de estabilidade corresponde ao anulamento 


TECNICA N: 374 


da área a tracejado, condição que é traduzida 
pelo sistema de equações 


5) vsinM-—sin %=0 


6) (—3o) sin do +» (cos dy — cos do) =0. 


Estas relações, para um dado valor de 2%, per- 
mitem determinar os valores de d; e », o que fiz 
por um processo iterativo, utilizando um cal- 
culador digital. 

Na fig. 2 traduz-se a relação entre o ângulo de 
potência inicial, do, e 31, ângulo de potência 
máximo durante a oscilação subsequente à per- 
turbação, e a tensão limite de estabilidade (tradu- 
zida pela relação v), para duração indefinida da 
perturbação de tensão. 


0 | T Ss tri) A 


À ms 0 
O 720 3% 4a SO 6 TO BO IO 106 igrau) 


Fig 2 — Relação entre o ângulo de potência inicial, 5,, 

e a tensão limite da estabilidade (traduzida pela razão v) 

podendo ser mantida indefinidamente, e relação entre 5, 

e o ângulo de potência máximo, 51, para essas condições 
limites. 


Para valores de » inferiores ao limite atrás in- 
dicado, a condição de estabilidade implica uma 
duração da perturbação de tensão inferior a deter- 
minado valor máximo. 

Como se verifica facilmente, a condição de 
estabilidade implica que seja positiva a área 
(algébrica) a tracejado na fig. 3, sendo 24 o ân- 
gulo de potência no instante em que se dá o 
restabelecimento da tensão. 

A duração máxima correspondente ao limite 
de estabilidade é traduzida pela condição 


(=— 29) sin 2,— (1 +) cos 9, 


7) cos 4= e 
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Fig. 3 — Relação entre o ângulo de potência e a potência 

trocada com a rede, P, antes e depois da perturbação de 

tensão (»==1),e durante a perturbação de tensão (» =="). 

Representação correspondente à verificação de estabilidade 

para uma perturbação de duração finita com instabilidade 
no caso de duração infinita. 


Tem-se assim, sob forma explícita, ?, em fun- 
ção de 9, e de v. 

A determinação da duração correspondente a 
ser atingido o ângulo de potência à, implica a 
integração da equação diferencial do movimento, 
designadamente traduzida pelas relações 1) e 3). 

Para dar maior generalidade à análise numé- 
rica, tem interesse uma mudança da variável 
tempo por uma variável auxiliar 7, tal que 


8) d+ = Vad = 


Esta mudança de variável corresponde no 
fundo a uma mudança de escala de tempos, 
tomando como unidade de tempo 
9) 0 = AT, 

Nesta unidade, a equação diferencial do movi- 
mento tem como únicos parâmetros 9, e », rela- 
cionando, para um par de valores desses pará- 
metros, ) e 7 

A unidade de tempo “ pode ser interpretada 
como a média geométrica entre o tempo que no 
movimento síncrono corresponde a um desloca- 


191 


mento eléctrico de um ângulo 2 radiano e o 
tempo igual à relação entre a energia cinética, W, 
associada ao rotor da máquina, e a potência 
«máxima», P, , que ela pode trocar com a rede, 
para f.e.m. equivalente E, tensão V, e reactân- 
cia equivalente X. 

No caso de uma frequência 50 Hz tem-se 


10) 0 = 0,0798 5! 


. V Po , 


e no caso de uma frequência 60 Hz têm-se 
: e /W 
11) 1=0,0728 8 4/5. 
E 


No caso de a energia cinética dominante ser 
a do rotor da máquina síncrona própriamente 
dita e para os valores correntes dos parâmetros, 
tem-se que O é da «ordem de grandeza» de 0,15. 

A equação diferencial do movimento foi inte- 
grada para valores de 9, de O a 90º, por interva- 
los de 5º, e para valores de v de 1 a O mas inferio- 
res ao valor limite de » correspondente a estabi- 
lidade para duração indefinida da diminuição de 
tensão. 

Na fig. 4 representam-se os resultados obtidos. 

É interessante salientar que para valores 7 


d 
| 


03 04 05 06 07 08 09 1 7 


oO 01 02 


Fig. 4 — Relação entre 3,,, v e 7 correspondente 
às condições limites de estabilidade. 
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elevados (superiores a cerca de 10 para 2, = 10º) 
as curvas de parâmetro 2, num diagrama (v, 7) 
coincidem sensivelmente com as respectivas 
assimptotas, isto é, para durações elevadas da 
perturbação de tensão e estabilidade é funda- 
mentalmente condicionada pela tensão durante 
a perturbação. 

Tem interesse caracterizar as perturbações de 
tensão por um parâmetro simples, que, aproxi- 
madamente, se possa associar à «gravidade» da 
perturbação. 

Considerei sistematicamente, para as condições 
limites de estabilidade, os seguintes parâmetros : 


5% 
12) Pa = (1 — v?) dr 
jN 
13) Pp = (1=v)º dr 
y 
14) p=) (1-9) ds 
IA 
AT 
15) Pa = dr 
(1—»v) dr 
16) Pe ei - TO 
dr 


Estes parâmetros foram calculados na escala 
de tempos 7, caracterizada por uma unidade de 
tempo “. Para a sua conversão a outra unidade 
de tempo, basta notar que os valores numéricos 
dos quatro primeiros são inversamente propor- 
cionais à unidade de tempo e o do último é 
independente da unidade de tempo. Designo os 
valores correspondentes, mas utilizando em vez 
de dz, o tempo dt, tomando como unidade o 
segundo, respectivamente por P,, P,, P., Py P.. 

Nas figs. 5,6,7,8e 9 representam-se estes 
parâmetros para condições limites de estabilidade, 
em função de 3, e ». 

Sob o ponto de vista restrito de existência 
ou não de estabilidade e para o caso de uma 
mesma máquina, caracterizada por um só valor 
de O, as condições limites correspondentes a 
um mesmo valor de 3, traduzem «igual» gravi- 
dade da perturbação. Como se verifica fácil- 
mente, nenhum dos parâmetros p,, Pp, Pc, Pq 
e p. satisfaz exactamente à condição de «unici- 
dade» para «igual» gravidade. 
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No entanto, como é óbvio, convém que, se 
possível, o parâmetro que «caracterize» a gra- 
vidade da perturbação de tensão seja «homoté- 
tico» para máquinas diferentes, eventualmente 


Fig. 5 — Relação entre 3,, v e o parâmetro p, correspon- 


dente à perturbação de tensão, para condições limites de 
estabilidade. 
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Fig. 6 — Relação entre do“ e o parâmetro Pp correspon- 


dente à perturbação de tensão, para condições limites de 
estabilidade. 


TÉCNIÇA N.º 374 


caracterizadas por valores de ? distintos. Desta 
forma, dois índices «iguais» sob o ponto de 
vista de uma máquina, caracterizada por um 
valor de 9 e outro de à, , deverão ser iguais sob 
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Fig. 7 — Relação entre 3, » e o parâmetro p, correspon- 


dente à perturbação de tensão, para condições limites de 
estabilidade. 
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Fig. 8 — Relação entre d,,v e o parâmetro py correspon- 


dente à perturbação de tensão, para condições limites de 
estabilidade. 
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Fig. 9 — Relação entre 3,, e o parâmetro Pe correspon- 


dente à perturbação de tensão, para condições limites de 
estabilidade. 


o ponto de vista de outra máquina, caracteri- 
zada por valores diferentes de 9 e de 5,. Este 
condicionamento desaconselha o uso do parâme- 
tro à, correspondente a limite de estabilidade, 
dados 7 e 9. 

Convém naturalmente, por outro lado, que 
um parâmetro que caracterize a perturbação 
tenha forma simples e por outro que, se possi- 
sivel, seja linear em 7. 

Da análise das figs. 5 a 9 conclui-se que, 
embora de maneira imperfeita, os parâmetros 
Pa € Pp; podem caracterizar satisfatóriamente o 


efeito da perturbação de tensão sob o ponto de 
vista de estabilidade dinâmica, para pequenas 
durações da perturbação. Para perturbações de 
longa duração, em que, como vimos atrás, é 
dominante a amplitude da perturbação, é mais 
significativo o parâmetro p.. 

O parâmetro p, é proporcional à diminuição 
provocada por uma perturbação de tensão na 
energia consumida numa carga caracterizada por 
uma impedância constante, portanto com elas- 
ticidade 2 em relação à tensão. O parâmetro p, 
é proporcional à diminuição de tensão numa 
carga caracterizada por uma elasticidade 1 em 
relação à tensão. 

Para concretização, sob o ponto de vista glo- 
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bal de um conjunto de máquinas de uma rede, 
do critério de escolha dos parâmetros caracte- 
rísticos da perturbação, suponhamos que os 
valores de à, dessas máquinas têm os valores 
5º; 10º; 25": 20"; 25 ; 30º; 35º; 40"; 45; 50”, 
e que os valores de O têm os valores 0,05.a; 
0,1.a; 0,15.a; 0,2.a; 0,25.a, sendo igual a pro- 
porção de máquinas para cada uma das diferen- 
tes 50 combinações de 9, e 1. 

Seja n a proporção de máquinas deste con- 
junto que perde o sincronismo com uma pertur- 
bação caracterizada pelos valores T (duração) 
e v (amplitude relativa) durante a perturbação. 

Nas fig. 10, 11 representa-se respectivamente 
a relação entre P, en e entre P. e 1. 


+. 


7 + 1p 2 
Fig. 10 — Relação entre o parâmetro P, e 1, proporção 
das máquinas que perdem o sincronismo, num exemplo 
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Fig. 11 — Relação entre o parâmetro P. e «, proporção 
das máquinas que perdem o sincronismo, num exemplo. 
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Da análise destas figuras conclui-se que, para 
perturbações de pequena duração, há uma corres- 
pondência aproximadamente biunívoca entre P. 


e n, e para perturbações de longa duração entre 
P. e, isto é, P. é adequado para «índice» da 
gravidade de perturbações de pequena duração e 
P. para índice da gravidade de perturbações de 


«longa» duração. A «transição» entre «pequena» 
e «longa» duração, para este efeito, dá-se para 
durações da ordem 10.4, ou seja da ordem de 1s 
para os valores «usuais» dos parâmetros das 
máquinas. 

Em conclusão, quanto a estabilidade de máqui- 
nas síncronas dentro de um ponto de vista glo- 
bal e simplista para o conjunto de uma rede, po- 
dem considerar-se como parâmetros adequados 
para traduzir a gravidade de uma perturbação de 
tensão: 


a) Para perturbações de curta duração (inferior 
a cerca de 15) 


ê 
17) P, = ) (1—v) dt, 
sendo T a duração da perturbação ; 


b) Para perturbações de longa duração (supe- 
rior a cerca de 15) 


1 T 


Com um pouco mais de generalidade, pode 
dizer-se que a zona de separação entre as dura- 
ções «curtas» e «longas» atrás indicadas é da 
ordem de 5 a 10 vezes o valor «médio» de 9 nas 
máquinas em causa. 


4. RELAÇÕES ADOPTADAS PARA O PRO- 
BLEMA DA SOBREINTENSIDADE ASSO- 
CIADA À PERTURBAÇÃO DE TENSÃO 


Dentro das hipóteses indicadas em 1), tem-se 
que: 


a) Durante a perturbação de tensão é (!) 
19) X P=E+ V— 2EVcos 5. 
b) Após a perturbação de tensão é (!) 


20) XP=E4+VI-2EV,cosi, 
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c) Na hipótese de não haver perturbação é (!) 
sa 2 2 
21) X L=E + V,—2E V, cos 9,. 


Como se verifica facilmente, supondo uma 
duração suficientemente curta para que se possa 
desprezar a transmissão de calor dos enrolamen- 
tos do induzido ao exterior, o aumento de tem- 
peratura de um enrolamento é apenas função de 


to 
22) G=| Fdt. 
ti 


Com efeito, nessas hipóteses, tem-se que 


dn En) TF 
nai dt g(n) 
sendo n a temperatura, t (1) a resistência (!) e g (1) 
a capacidade calorífica (!) de um troço de condutor, 


donde 


24) q (n,t,t) =[ "80 qu=[ “ra. 
tu E(n) tu 


Portanto, para um par de valores (t,, t.), a 
função 7(n) define o valor da temperatura, e é 
por sua vez igual a G. 

Com um certo erro, pequeno para fortes sobre- 
intensidades de curta duração, pode ainda admi- 
tir-se que o aumento de temperatura correspon- 
dente a essa sobreintensidade I, e em relação à 
temperatura estabilizada com uma corrente I,, é 
função de 


t 
25) H=[ (BI!) dt. 
t 


Considerando para origem dos tempos o iní- 
cio da perturbação de tensão, seja 


? t 
dd Qdo, É [ema 


9 “EV o 
X ota 
o E =) ds. 
o o 
Vem: 
a) parat< T 


27) t= f [$e- 1)+2 (cosê, — vcos » | ds; 


(') Consideram-se valores de corrente referidos a um 
esquema monofásico equivalente. 
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b) para t>T 


aa 
28) € =[ | < (2-1) + 2 (cos? ,-vcos ) [d+ 


t 
+f 2 (cos 3,— cos 3) dr. 
T 


dd 


A grandeza É pode ser interpretada como a 
relação entre o acréscimo de energia reactiva, 
em relação ao valor anterior à perturbação, e o 
produto da constante de tempo Ú pela potência 
máxima (para E, V, e X constantes e 9 arbitrá- 


rio) que a máquina síncrona pode trocar com a 
rede. 

Esta grandeza depende do intervalo de tempo 
a que se refere, e das condições anteriores à per- 
turbação através dos parâmetros : e 9,. 


5. ANÁLISE NUMÉRICA DA SOBREINTEN- 
SIDADE NO INDUZIDO 


Para efeitos de aquecimento, interessa desig- 
nadamente conhecer o valor de É (ou de gran- 
deza similar), durante a perturbação de tensão e 
no período subsequente em que a corrente se 
mantém elevada. 

Para uma análise sistemática elaborei um pro- 
grama que, com base na integração numérica 
das equações 1 e 3,e para os diversos valores 
das grandezas, calcula os valores de É correspon- 
dentes aos seguintes intervalos de tempo: 


a) Durante a perturbação de tensão (=). 


b) Desde o início da perturbação até ao ins- 
tante T,, posterior à perturbação de ten- 
são, em que a corrente atinge novamente 
o valor anterior à perturbação ou é infe- 
rior a este último. Após T, a corrente passa 
a ser inferior a I,, seguindo-se uma osci- 
lação de I em torno do valor 1, (=). 


O critério de escolha do intervalo T, para cál- 
culo de É, é naturalmente um tanto arbitrário. 
Afigura-se no entanto simples e representativo 
para uma análise aproximada dos efeitos de uma 
perturbação de tensão sob o ponto de vista de 
sobreintensidade no induzido de máquinas sín- 
cronas, dentro das hipóteses feitas. 

Nas figs. 12, 13 e 14 apresenta-se um extracto 
dos resultados obtidos, indicando-se, para E = 
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= 1,l,“,e T, em função de T, duração da per- 
turbação de tensão, para 9, = 10º; 20º; 30; 40"; 
50º. A análise feita engloba apenas condições 
correspondentes à existência de estabilidade diná- 


mica das máquinas. 


Fig. 12 — Relação entre a duração, T, da diminuição de 

tensão, de amplitude relativa », e o parâmetro Z,, repre- 

sentativo do integral do acréscimo do quadrado de cor- 
rente durante a diminuição de tensão. 
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Da análise destas figuras, conclui-se que, sob 
um ponto de vista global e de primeira aproxi- 
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Fig. 12 (cont.) —Relação entre a duração T, da diminuição 

de tensão, de amplitude relativa v, e o parâmetro 1,, re- 

presentativo do integral do acréscimo do quadrado de 
corrente durante a diminuição de tensão. 


c) 36 = 30º 
a) 3 = 10º 
Fig. 13 — Relação entre a duração T, da diminuição de Fig. 13 (cont.)-Relação entre a duraçãoT, da diminuição de 
tensão, de amplitude relativa », e o parâmetro 1, repre- tensão, de amplitude relativa v, e o parâmetro %,, repre- 
sentativo do integral do acréscimo do quadrado de cor- sentativo do integral do acréscimo do quadrado de cor- 
rente até o ângulo de potência atingir novamente o valor rente até o ângulo de potência atingir novamente o valor 
inicial, 56» inicial, 3, 
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Fig. 13 (cont.) -Relação entre a duração T, da diminuição de 

tensão, de amplitude relativa v, e o parâmetro 1, repre- 

sentativo do integral do acréscimo do quadrado de cor- 

rente até o ângulo de potência atingir novamente o valor 
inicial, 56 » 


mação, o índice P. pode ser considerado «repre- 
sentativo» da gravidade de perturbação de ten- 
são quanto à sobreintensidade resultante no in- 
duzido de máquinas síncronas, na medida em 
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se seja legítimo considerar as hipóteses indicadas 
em 1). 
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b) do == 20º 


Fig. 14 — Relação entre a duração T, da diminuição de 
tensão, de amplitude relativa v, e o tempo, T,, de restabe- 
lecimento do ângulo de potência inicial, 5, . 
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d) 5, = 40º 


Fio. 14 (cont,)— Relação entre a duração T, da diminui- 
ção de tensão, de amplitude relativa v, e o tempo, T,, de 
restabelecimento do ângulo de potência inicial, 3, « 
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e) Jo = 50º 


Fig. 14 (cont.) — Relação entre a duração T, da diminui- 
ção de tensão, de amplitude relativa v, e o tempo, T,, de 
restabelecimento do ângulo de potência inicial, 5, - 


Tem aqui lugar um esclarecimento do motivo 
da «discontinuidade» da relação 


T=T(T. 


Na fig. 15 representa-se a evolução de à em fun- 
ção do tempo para diversos valores da duração 
da redução de tensão. Após a redução de tensão 
a máquina entra em movimento oscilatório em 
torno da frequência de sincronismo, variando à 
aproximadamente entre d, e um valor máximo 
dm» com um «período» T.. 


Para 
E E Mo 


T;—T 
é função decrescente de T, sendo 


a diferença 


29) lim T;=T. 


T-—s T —b 
o 


Aumentado T, para 
T=T,+4o0 


dá-se uma discontinuidade no valor de T,, pois, 
ao restabelecer-se a tensão, inicia-se uma oscila- 
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Fig. 15 — Evolução do ângulo de potência, 5, em função do tempo, para vários valores 
de duração da diminuição de tensão, T. 


ção de d em fase de crescimento a partir de um 
valor superior a 2,, só sendo atingido o valor 2, 
após uma «oscilação». 

Discontinuidades análogas verificam-se, Obvia- 
mente, para 


30) [=kT,, 


com k inteiro positivo e não nulo. 

Claro que o amortecimento das oscilações, não 
tido em conta no modelo considerado, influi 
significativamemente no valor de à após um 
número elevado de períodos. 


6. CONCLUSÃO 


Em condições em que sejam aproximadamente 
válidas as hipóteses indicadas em 1), os resultados 
numéricos apresentados permitem averiguar as 
condições de estabilidade de máquinas síncronas 
perante uma perturbação de tensão e as conse- 
quências dessa perturbação sob o ponto de vista 
de sobreintensidade no induzido. 

Dentro da validade dessas hipóteses, quanto a 
estabilidade e sob um ponto de vista global do 
conjunto de máquinas síncronas de uma rede, 
tem-se que: 


a) Para perturbações de tensão de «curta» 
duração, a gravidade da perturbação é satis- 
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fatôriamente traduzida pelo parâmetro P., 
sendo 


+T 


31) P no (1 —»v) dt. 


c 


b) Para perturbações de «longa» duração, a 
gravidade da perturbação é aproximada- 


mente traduzida pelo parâmetro 


32) ARE ro 


vo 


(1=—v dt. 


Para a ordem de grandeza usual dos parâme- 
tros das máquinas e para este efeito, a transição, 
aliás não rígida, entre duração curta e longa 
verifica-se para tempos da ordem de 1 s. 

Semelhantemente, e para análogas hipóteses e 
ponto de vista e quanto a sobreintensidade no 
induzido de máquinas síncronas associadas à 
componente directa da tensão, o parâmetro P,. 
pode considerar-se um índice significativo da 
gravidade da perturbação de tensão. 

Afigura-se pois de interesse que os elementos 
estatísticos relativos a perturbações de tensão 
nas redes indiquem designadamente os parâme- 
tros P. e P,, ou qualquer par das três gran- 


dezas P., P. e T, dado que elas estão relacio- 
nadas por 


33) P.=P,.T. 
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APLICAÇÃO DOS COMPUTADORES AO CÁLCULO 
DINÂMICO DE ESTRUTURAS Ú 


RESUMO 


Situa-se a análise dinâmica das estruturas relativa- 
mente aos vários domínios a que está ligada, 

Enumeram-se os seus objectivos e referem-se as caracte- 
résticas formais dos métodos de que se dispõe para os 
atingir. 

Realça-se a contribuição prestada pelos computadores 
etectrônicos digitais em todos os campos da análise dinã- 
mica, tanto ao nível de investigação como de projecto, 


1 — INTRODUÇÃO 


O estudo do comportamento de estruturas 
actuadas por solicitações em regime dinâmico 
tem sofrido um desenvolvimento considerável 
nos últimos vinte anos, constituindo presente- 
mente uma disciplina perfeitamente caracterizada 
dentro da Teoria das Estruturas. 

Dentre as causas a que se deve este desenvol- 
vimento, destacam-se as seguintes: 


— O progresso dos conhecimentos nos vários 
campos ligados à análise dinâmica, designada- 
mente em Estatística, Teoria das Vibrações, Sis- 
mologia, Aerodinâmica, Resistência dos Mate- 
riais e Teoria das Estruturas, permitiu reunir um 
conjunto de elementos necessários ao estudo dos 
problemas em causa ; 

— a evolução das estruturas torna indispensá- 
vel considerar no respectivo dimensionamento 
solicitações cujo carácter dinâmico assume grande 
importância ; é o caso de torres, chaminés, edi- 
fícios de grande altura e coberturas de grande 
vão, expostas à acção do vento e dos sismos; 

— a instalação de maquinaria em instalações 
industriais exige o dimensionamento das estru- 
turas sob a acção das solicitações transmitidas 
pelas máquinas ; 

—e finalmente, o recurso a meios de cálculo 
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muito potentes, ou seja, a computadores electró- 
nicos digitais permite presentemente a análise 
dinâmica de estruturas de grande complexidade. 

A análise dinâmica de uma estrutura com- 
preende as seguintes fases : 


— definição da solicitação dinâmica ; 

— definição de um modelo matemático que 
represente sob os pontos de vista de deformabi- 
lidade e absorção de energia o comportamento 
da estrutura ; 

— aplicação da teoria das vibrações ao estudo 
do comportamento do modelo matemático. 


O objectivo da análise dinâmica é a determi- 
nação dos deslocamentos, velocidades e acelera- 
ções adquiridos pelos elementos da estrutura. 
A partir da deformada dinâmica obtida, o cál- 
culo dos esforços prossegue nos moldes habi- 
tuais da Teoria das Estruturas. 


2 — SOLICITAÇÕES EM REGIME DINÂMICO 


O carácter estático ou dinâmico de uma solici- 
tação é relativo à estrutura sobre a qual esta 
actua. Assim, uma dada solicitação tem carácter 
dinâmico em relação a uma estrutura quando as 
forças de inércia que se desenvolvem nesta têm 
valores significativos em relação às restantes 
forças que intervêm no equilíbrio da estrutura. 

A mesma solicitação pode portanto actuar em 


(*) Comunicação apresentada às II Jornadas Luso-Brasileiras de Engenharia Civil — Rio de Janeiro e S. Paulo, 1967. 
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regime estático ou em regime dinâmico sobre 
estruturas de características diferentes. 

As solicitações dinâmicas podem classificar-se 
em determinísticas e aleatórias, conforme é 
conhecido o seu valor num dado instante ou 
apenas a respectiva distribuição estatística. As 
solicitações determinísticas dividem-se por sua 
vez em periódicas e aperiódicas [1]. 


2.1 — Solicitações deterministicas periódicas 


Pertencem a este grupo as solicitações origi- 
nadas pelo funcionamento em regime permanente 
de máquinas de vários tipos. A representação 
analítica das solicitações varia muito conforme 
o tipo de máquina. Assim, as máquinas rotativas 
de alta velocidade, tipo turbinas, alternadores, 
centrifugadores, motores eléctricos, caracteriza- 
dos por regimes de carga sensivelmente unifor- 
mes, dão origem a solicitações sinusoidais sim- 
ples; as máquinas alternativas e certos tipos de 
máquinas-ferramentas produzem solicitações si- 
nusoidais de várias frequências (fundamental e 
harmónicas); as máquinas de percussão, tipo 
martelos-pilões, saca-bocados, etc., dão origem a 
solicitações correspondentes a impulsos de curta 
duração. 


2.2 — Solicitações determinísticas aperiódicas 


Pertencem a este grupo as solicitações devidas 
ao regime transitório (arranque e paragem) de 
máquinas de qualquer dos tipos indicados bem 
como as explosões na atmosfera e a maioria das 
solicitações devidas ao tráfego de veículos. 

Relativamente às explosões interessa distinguir 
as explosões na atmosfera, em que a representa- 
ção analítica da onda de pressão é bem conhe- 
cida, das explosões subterrâneas, que dão origem 
a movimentos do solo de natureza muito com- 
plexa, semelhante à dos fenómenos sísmicos, a 
que se fará uma breve referência em seguida. 


2.3 — Solicitações aleatórias 


Estão neste grupo as solicitações devidas ao 
vento, sismos, certos tipos de máquinas e de 
tráfego. 

Caracteriza-se este tipo de solicitações por meio 
da respectiva distribuição estatística. Relativa- 
mente ao vento e aos sismos, esta distribuição 


202 


refere-se à velocidade das rajadas u (t) e às ace- 
lerações no terreno a” (t) e admite-se que possui 
as seguintes características : 


— é estacionária num intervalo de tempo con- 
venientemente escolhido ; 
—tem uma distribuição de amplitude gaus- 
seana ; 
— tem um valor médio nulo no caso dos sis- 
mos e diferente de zero no caso do vento; 
— a repartição de energia pelas diversas fre- 
quências, caracterizada pelo espectro de potência, 
tem a forma indicada na fig. 1, [1] sendo afectada 
pelas condições locais. 


É de notar que a influência das condições lo- 
cais sobre o espectro de potência é quase desco- 
nhecida no que respeita aos sismos, ao contrá- 
rio do que sucede para o vento. 

Na figura 1 indicam-se também valores carac- 
terísticos da frequência própria de vários tipos 
de estruturas, o que permite avaliar a importân- 
cia do carácter dinâmico do vento e dos sismos 
sobre cada um dos tipos indicados. Assim, veri- 
fica-se que o vento só constitui uma solicitação 
de carácter dinâmico para estruturas de baixa 
frequência própria, enquanto os sismos atingem 
especialmente estruturas de frequências próprias 
mais elevadas. Fica deste modo evidenciada a re- 
latividade do carácter estático ou dinâmico de 
uma solicitação. 


3 — GENERALIDADES SOBRE EQUILÍBRIO 
DINÂMICO 


A condição de equilíbrio de uma estrutura su- 
jeita a solicitações em regime estático e dinâmico 
é o anulamento da resultante das solicitações 
actuantes em cada um dos seus elementos. 

No caso de todas as solicitações aplicadas se- 
rem forças, deverá verificar-se, em cada elemento 
da estrutura, equilíbrio entre todas as forças 
actuantes, que podem ser dos seguintes tipos: 


F. — forças de inércia 

F. — forças dissipativas de atrito ou amorteci- 
mento 

F. — forças de restituição ou ligação 

F. — forças aplicadas 


ou seja, deverá ser em cada elemento da estru- 
tura : 
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DENSIDADE ESPECTRAL REDUZIDA, e e S (1) 
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PERIODO (5) 
Fig. 1 — Carácter dinâmico do vento €« dos sismos 
ELERTHE+rHE=Oo (1) | o estudo do comportamento de osciladores linea- 


A equação (1) exprime o teorema de d'Alem- 
Lert. 

Os termos F e F, dependem da aceleração e 
da velocidade adquiridas pelos elementos da 
estrutura; se forem desprezáveis, a equação (1) 
reduz-se a: 


EF + F,=0 (2) 


obtendo-se uma condição de equilíbrio estático. 

O número de equações que descreve o movi- 
mento de uma estrutura é o respectivo número 
de graus de liberdade. Uma estrutura contínua 
tem infinitos graus de liberdade; é usual discre- 
tizá-la, considerando a massa concentrada num 
número finito de pontos. Considera-se as forças 
aplicadas actuando únicamente nesses pontos, 
obtendo-se assim um sistema com um número 
finito de graus de liberdade. 

O movimento oscilatório ou vibração é carac- 
terizado por as forças de restituição contrariarem 
constantemente o movimento, tendendo a con- 
duzir a estrutura a uma posição de repouso. 
Estas forças podem ou não ser proporcionais 
aos deslocamentos, dando origem a oscilações 
lineares ou não lineares. 

O objectivo da teoria das vibrações é portanto 
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res e não lineares, com um número finito ou infi- 
nito de graus de liberdade, sob a acção de soli- 
citações dinâmicas de qualquer dos tipos indi- 
cados em 2. 


4 — MÉTODOS GERAIS DE ANÁLISE 
DINÂMICA 


4.1 — Oscilador linear de 1 grau de liberdade 


Assim, em qualquer dos casos da fig. 2, 
pode supor-se a massa da estrutura, m, concen- 
trada num ponto e estudar-se as oscilações de 
translacção que serão regidas por uma única 
equação de equilíbrio. Admitindo que o amorte- 
cimento é do tipo viscoso, ou seja, proporcional 
à velocidade, os termos da equação (1) tomam 
a forma: 


F;= cx 
F. = Kx — em regime linear 


F.== F(t) no caso (a) 
F.=— m a (t) nos casos (b) e (c) 


e a equação diferencial do movimento é 
mx +cx+ Kx= E (3) 
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Máquina 


a(t) 


a)Transiacção vertical de 
fundação de máquina 


b) Oscilação sísmica de 
reservatório elevado 


Cunjunto muito 
rigido 


Conjunto deformávet 


C) Oscilação sísmica de 
estrutura de edificio 


Fig. 2 — Osciladores de 1 grau de liberdade 


sendo : 


c — coeficiente de amortecimento viscoso 

K — coeficiente de rigidez ou simplesmente 
rigidez 

x — elongação (absoluta ou relativa ao movi- 
mento do solo, conforme se tratar dos 
casos (a) ou (b) e (c)). 


A integração da equação (3) permite, como é 
bem sabido, determinar a frequência própria 
mm, = (K/m)!/» e a resposta do oscilador, x (t), 
em regime livre (F, = 0) e forçado (F. & 0), 
sujeito à actuação de solicitações determinísticas 
e aleatórias de qualquer dos tipos referidos em 
3. fas. 

Os esforços que se desenvolvem na estrutura 
devido à oscilação são em seguida determinados 
considerando aplicadas as solicitações F, = Kx, 

Na definição da rigidez K é necessário ter em 
conta as eventuais alterações das características 
elásticas dos materiais devidas ao regime diná- 
mico. 

Em relação a problemas em que intervenha a 
oscilação do solo a rigidez tem normalmente que 
ser determinada por via experimental [4,5]. 


4 2 — Oscilador linear de vários graus 
de liberdade 


Na grande maioria dos casos, o modelo mate- 
mático mais adequado para representar o com- 
portamento dinâmico de uma estrutura é o osci- 
lador de vários graus de liberdade. 

O critério de localização e quantificação das 
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massas concentradas representa um problema 
delicado sempre que na estrutura não há de 
facto pontos de concentração de massas. No 
caso de análise dinâmica de estruturas de edifí- 
cios é corrente considerar a massa concentrada 
ao nível dos pisos. 

De um modo geral, o número de massas con- 
centradas que é necessário considerar para repre- 
sentar o sistema real é aproximadamente duplo 
do número de modos de vibração que contri- 
buem de forma significativa para a resposta do 
oscilador. 

A formulação das equações de equilíbrio no 
caso de osciladores lineares de vários graus de 
liberdade conduz a sistemas de equações dife- 
renciais ordinárias de 2.º ordem. A integração 
destes sistemas faz-se directamente ou decom- 
pondo-os em equações independentes por con- 
veniente mudança de coordenadas. Em qualquer 
dos casos, a escrita das equações é muito sim- 
plificada recorrendo à representação matricial. 
E assim, na prática corrente, os problemas de 
análise dinâmica de estruturas com um número 
finito de graus de liberdade são analisados por 
um dos três métodos seguintes : 


4.2.1 — Integração directa de sistemas de equações 
diferenciais 


Em problemas formulados através de um 
pequeno número de equações torna-se simples a 
respectiva integração directa. É este por exemplo 
o caso de oscilações lineares e de translacção do 
bloco de fundação, suposto assente sobre um 
meio elástico, indicado na fig. 3. 
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O movimento do bloco é constituído por uma 
translacção x(t) e por uma rotação O(t). O 
modelo matemático que o representa é portanto 


F(t) 


e——— > 


Fig. 3 — Oscilador de dois graus de liberdade 


um oscilador com dois graus de liberdade, cujas 
equações são da forma 


mxHEÉ +E=E 
, a Fr (4) 
I9 +N=0 


sendo a segunda equação relativa à oscilação 
angular. 

A integração do sistema (4) permite obter as 
frequências próprias e modos de vibração do 
oscilador, bem como a resposta no tempo, x (t) 
e 9(t). Em problemas deste tipo as solicitações 
aplicadas são geralmente determinísticas, de 
forma simples, pelo que o número de cálculos a 
efectuar é reduzido. 


4.2.2 — Análise modal 


+ 


Este método é sem dúvida o mais potente e 
versátil em análise dinâmica pois permite o tra- 
tamento sistemático de qualquer tipo de oscila- 
dor linear sujeito a qualquer tipo de solicitação. 
Em linhas gerais consiste no seguinte [6, 7]: 


Seja a estrutura indicada na fig. 4, sujeita à 
actuação de forças horizontais em regime diná- 
mico. Considera-se a massa concentrada ao nível 
dos pisos e pretende estudar-se as oscilações de 
translacção horizontal x (t'. 

A estrutura é redutível ao oscilador de 3 graus 
de liberdade indicado na fig. 4. 

As equações de equilíbrio, dada a linearidade 
do sistema, podem escrever-se sob a forma 


m, X, on C'x, + SK, x= | Lj=1,2,3 (5) 
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ou, adoptando a representação matricial : 


[m] (x) + [C] (x) + [K] bi) = (E) (9) 
As eg. (6) são válidas para osciladores de n 
graus de liberdade, sendo: 


[m] — matriz das massas, diagonal, de ordem n. 

[C] —matriz de amortecimento, suposto de 
tipo viscoso. 

[K] — matriz de rigidez 


As frequências » e os modos de vibração (6) 


do oscilador correspondem aos valores próprios 


m3 


m1 


Fig. 4 — Oscilador de 3 graus de liberdade 


e vectores próprios da matriz [m”! K], ou seja: 


mw? — valores próprios de [m”! K] dados por 
| K -— nm? m | = 8) 
(6! == vectores próprios de [m-! K] 


É corrente calcular wo e [9] a partir não da 
matriz rigidez [K] mas da matriz de deformabi- 
lidade [4] sendo 

[4] = [KIM 
e: 
1/ == valores próprios de [4 m] 
(4) = vectores próprios de [4 m] 


Para o estudo da oscilação em regime forçado 
adopta-se a transformação linear de coordena- 
das x — : definida por 


Dx) = [9] (8) 


e dada a ortogonalidade dos modos de vibração 
(vectores próprios de matriz simétrica real) o sis- 
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tema de equação (6) transforma-se em n equações 
independentes, da forma 


| a TE 

+Hanve ta eo j=1,2,..n 

i E KA i (4) [m] EA (7) 
Cada equação refere-se a um modo de vibra- 

ção. n, designa o amortecimento relativo do 

modo i. Admite-se que [C] é diagonalizada pela 

transformação que diagonaliza [K], para o que 


deve ser [8]: 
[C]=a [m] + E [K] 


Designa-se por massa generalizada M, e por 
factor de participação R, relativos ao modo i 


os valores : 
M,= (6)" [ml (6; 
R, are [9,)" (E) 


A integração das equações (7) conduz aos va- 
lores de :; (t), a partir dos quais se calcula suces- 


sivamente: 
bj = [9] (9) 
[E = [KI (9 


obtendo-se os deslocamentos e as forças de res- 
tituição. A partir dos respectivos valores máximos 
calculam-se os esforços na estrutura. 

A determinação das séries cronológicas : (t), 
x (t), F(t) representa um grande volume de cál- 
culo e dos resultados obtidos só se utilizam os 
valores máximos, visto serem os que conduzem 
a maiores esforços na estrutura. 

Como por outro lado nem todos os modos con- 
tribuem de forma significativa para o movimento, 
é usual determinar apenas os valores máximos 
das amplitudes e forças em cada modo, sobre- 
pondo os resultados através de critérios de base 
estatística. Este processo é particularmente empre- 
gado para solicitações aleatórias, especialmente 
em Engenharia Sísmica |9). 

Se as solicitações aplicadas F (t) forem sinu- 
soidais de frequência w e se puder desprezar a 
influência do amortecimento, as eq. (6) tomam 
a forma 


[m] (x) + [K] |xb = (Ff sen ot (9) 


e sendo 


(8) 


x) = (A? sen » t 


obtém-se fácilmente 


JAP Ko? m] ! E) 
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Este caso é muito corrente em estruturas cujas 
fundações são actuadas por um regime de acele- 


rações sinusoidais a sen w t, devidas ao funcio- 
namento de máquinas instaladas nas proximida- 
des. É então 

J = e— a tn) 

Ez, a ms sen 0 t 


4.2.3. — Matrizes de transferência 


O conceito de matriz de transferência baseia-se 
nos trabalhos de Holzer [10], Myklestad e Thom- 
son [11 e 12] sobre análise de vibrações em vi- 
gas sujeitas a torção e flexão. Em linhas gerais, 
o método de Holzer consiste no seguinte: 

Seja um troço de veio cilíndrico com massas 
concentradas em vários pontos, conforme se in- 
dica na fig. 5. 

A rigidez de cada tramo de veio compreendido 
entre duas massas consecutivas é constante, po- 
dendo variar de tramo para tramo. 

Se o veio estiver sujeito a oscilações de tor- 
ção de frequência », a condição de equilíbrio da 
massa i será 

MM, 4 qui SM, E Lô, mt = 0 (10) 
em que: 


X — momento de torção 
I — momento de inércia polar da massa i 
9, — deslocamento angular da massa i 


Axi AXic1 


Xi 


ig. 5 — Veio cilíndrico não uniforme 


Sendo K a rigidez de torção em cada troço 


de veio será também: 
né 


MM. 
E o RE ES 
l Í 


1+t (11) 


i 


As equações (10) e (11) podem escrever-se 
sob a forma: 


1? 1 — (ud - 
EPE — IM, fe) I, O, 


Ry/8 12 
O, | == (3, + pass Es LM ( ) 
K 
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Se o veio tiver n massas concentradas, cons- 
tituindo um oscilador de n graus de liberdade, 
a partir das respectivas condições de ligação ao 
exterior, ou seja, dos valores de 9%, O, WN, e O, 
é possível, por iterações sucessivas, determinar 
os valores de w e 0),, ou seja, as frequências 
próprias e modos de vibração. O método é tam- 
bém aplicável ao estudo de oscilações forçadas 
sinusoidais desde que as solicitações aplicadas 
intervenham nas equações (12). 

As equações (12) podem ser escritas na forma 
matricial : 


di 2 I, 1 
(1 — ) needs | (E) 
K, K, - 


— 3 1, 1 


ou abreviadamente 


Zig = Eu), tl, (14) 

As matrizes coluna [z] chamam-se vectores de 
estado porque descrevem os esforços e as defor- 
mações em cada elemento da estrutura. [u], é a 
matriz de transferência do elemento i e os seus 
elementos dependem das características da estru- 
tura e da frequência de vibração. 

Os vectores de estado em dois elementos quais- 
quer de uma estrutura podem ser relacionados 
através das matrizes de transferência intermé- 
dias, visto ser: 


1Z), e [u], Zi, 
(z), = [u], (Z), = [u], [u), 42), 
(Z), = [0], [u],.,... [u], [u)] (2), = [PI AZ), 


Os elementos de |Z| e [u] dependem do pro- 
blema em causa. 

Myklestad e Thompson [11 e 12] desenvolve- 
ram um método análogo ao de Holzer para estu- 
dar a flexão, do qual deriva o conceito de matriz 
de transferência aplicado à flexão. 

O desenvolvimento do método das matrizes 
de transferência e a sua aplicação a um número 
muito variado de problemas deve-se a Pestel e 
Leckie [13], que tabelaram as matrizes de trans- 
ferência para um grande número de estruturas, 
o que permite a aplicação expedita do método 
na determinação de frequências próprias e modo 
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de vibração e em problemas de oscilação forçada 
devida a solicitações sinusoidais. 


4.3 — Oscilador linear com infinitos graus de 
liberdade 


O grau de aproximação dos métodos expostos 
depende do número de pontos em que se supõe 
concentrada a massa da estrutura. Se esta na 
realidade não apresentar pontos de concentração 
de massa, a discretização pode conduzir a um 
sistema cuja análise exija um volume de cálculo 
excessivo. E assim, torna-se em certos casos mais 
prático abordar o problema considerando o osci- 
lador com uma distribuição contínua de massa e 
elasticidade, ou seja, com um número infinito 
de graus de liberdade. 

As condições de equilíbrio são expressas neste 
caso por equações diferenciais de ordem igual ou 
superior a 2 que integradas definem as caracte- 
rísticas dinâmicas e a resposta da estrutura. 

A título de exemplo, considere-se o problema 
da oscilação de flexão de uma viga conforme se 
indica na fig. 6: 


F 


Fig. 6 - Vibração de uma viga sujeita a flexão 


Sendo 
2 3 4 
dx dx* dx* 
ar 
dx 


em que F, é a intensidade da força de restituição, 
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a condição de equilíbrio dinâmico será, despre- 
zando o amortecimento : 


em que: 


A — área da secção transversal 
Y — massa específica 


Obtém-se assim uma equação diferencial de 
4.º ordem cuja integração, entrando com as con- 
dições aos limites, relativas às ligações da viga 
ao exterior, permite obter as frequências próprias, 
modos de vibração e a resposta da estrutura, 
y (x, t). 

Se a solicitação aplicada ou a distribuição de 
massa e elasticidade apresentam descontinuidades 
é necessário introduzir condições aos limites adi- 
cionais que tornam o cálculo muito trabalhoso. 
Por este motivo os métodos baseados na teoria 
das vibrações de sistemas contínuos são normal- 
mente utilizados na análise dinâmica de estru- 
turas em que não há variações bruscas de carac- 
terísticas. 

Dentro destes métodos merece referência espe- 
cial o método desenvolvido por Kolousek [14] 
para a análise de estruturas recticuladas e arcos, 
sob acção de solicitações sinusoidais ou decom- 
poníveis nestas. , 

O método de Kolousek é formalmente análogo 
ao método dos deslocamentos utilizado em aná- 
lise estática. Assim, Kolousek deduziu e tabelou 
as fórmulas que dão os elementos da matriz rigi- 
dez de barras sujeitas a deslocamentos sinusoidais 
aplicados nas extremidades; procedeu de igual 
modo em relação aos esforços desenvolvidos nas 
extremidades de barras sujeitas a solicitações 
sinusoidais concentradas e distribuídas. A partir 
destes elementos o cálculo da deformada em 
regime dinâmico e dos esforços correspondentes 
processa-se tal como no método dos desloca- 
mentos. 

O método de Kolousek permite considerar 
amortecimento de tipo viscoso, recorrendo neste 
caso a funções de variável complexa, técnica que 
aliás também é adoptada no método das matrizes 
de transferência. 

É de notar finalmente que a aplicação dométodo 
dos deslocamentos em análise dinâmica pode ser 
feita através de processos iterativos [15], aná- 
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logos aos utilizados nos métodos de Cross e de 
Kani. 


4.4 — Comportamento em regime não linear 


É no campo da Engenharia Sísmica que a aná- 
lise dinâmica de estruturas em regime não linear 
apresenta maior interesse prático, visto que um 
sismo de forte intensidade produz nas estruturas 
correntes de edifícios oscilações de amplitude 
superior aos máximos admissíveis em regime 
elástico. 

As características de não linearidade de uma 
estrutura são expressas através do respectivo 
diagrama forças-deslocamentos. Indicam-se na 
fig. 7 alguns dos diagramas adoptados para es- 
truturas metálicas e de betão armado. Os dois 
primeiros são os mais correntes. 

As condições de equilíbrio são agora expressas 
por sistemas de equações diferenciais não linea- 
res e o método geralmente empregado para a sua 
resolução é a integração «passo a passo”, que em 
linhas gerais consiste no seguinte : 


— decompor a solicitação aplicada em impulsos 
de pequena duração ; 

— considerar que durante a actuação de cada 
impulso a estrutura se comporta em regime 
linear ; 

— a partir das condições iniciais do movimento 
calcular o acréscimo de amplitude, velocidade 
e aceleração nos elementos da estrutura de- 
vidos ao primeiro impulso ; 

— adicionar os valores obtidos aos iniciais ; 

— considerar a característica de restituição cor- 
respondente aos valores obtidos ; 

— calcular novamente os acréscimos de ampli- 
tude, velocidade e aceleração devidos ao im- 
pulso seguinte ; 

— proceder assim sucessivamente para cons- 
truir as séries cronológicas dos valores da 
amplitude, velocidade e aceleração que des- 
crevem o comportamento dinâmico da es- 
trutura. 


A formulação matricial continua a estar indi- 
cada neste tipo de problemas, desde que se con- 
siderem as alterações da matriz rigidez devidas 
ao comportamento não linear [18, 19). 

De entre os métodos específicos para o estudo 
de osciladores não lineares destacam-se dois : 


TÉCNICA N.º 374 


a) Diagrama elasto-plástico 


- F 


C) Diagrama de Jennings [16] 


b) Diagrama bi-linear 


d) Diagrama elasto-plástico de rigidez 
decrescente 117 


Fig. 7 — Diagramas forças-deslocamentos em regime não-linear 


— método gráfico do plano de fase |2], aplicá- 
vel a osciladores de um grau de liberdade 
sob a acção de solicitações determinísticas, 
especialmente indicado para o estudo de os- 
cilações em regime livre; 

— método de Ritz [20 analítico, baseado no 
cálculo variacional, aplicável a osciladores 
não lineares de um ou vários graus de li- 
berdade |21!, em regime forçado. 


5 — APLICAÇÃO DOS COMPUTADORES AO 
CÁLCULO DINÂMICO 


A aplicação numérica dos métodos da teoria das 
vibrações a estruturas complexas representa um 
volume de cálculo que só o advento dos compu- 
tadores electrónicos tornou praticável. Assim se 
compreende que o desenvolvimento da análise 
dinâmica de estruturas esteja intimamente asso- 
ciado ao desenvolvimento dos computadores e 
que se tenham explorado principalmente os mé- 
todos de formulação simples que envolvam grande 
número de operações de carácter repetitivo. 

Deste modo, na análise dinâmica de estruturas 
em regime linear adopta-se geralmente o método 
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da análise modal para osciladores comn graus de 
liberdade, o que implica a realização dos seguin- 
tes tipos de cálculo: 


— deformações devidas a forças unitárias, para 
determinação da matriz de deformabilidade 
[4]; 

— inversão da matriz de deformabilidade, para 
cálculo da matriz rigidez, [K]; 

— valores próprios e vectores próprios de 
[m-! K] ou [4 m], sendo [m] a matriz das 
massas, supostas concentradas em n pontos, 
para determinação das frequências próprias 
e modos de vibração ; 

— factores de participação e massa generali- 
zada em cada modo; 

— integração de equações diferenciais lineares 
de 2.2 ordem; 

— operações sobre matrizes e vectores. 


Todas as operações descritas correspondem a 
problemas correntes, para os quais existem pro- 
gramas nas bibliotecas do Centro de Cálculo, 
excepto a determinação da matriz [4] que cons- 
titui um problema básico em Teoria de Estruturas. 
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Pode assim afirmar-se que é actualmente pos- 
sivel efectuar a análise dinâmica de qualquer 
estrutura desde que se saiba determinar a res- 
pectiva matriz de deformabilidade, ou seja, desde 
que se saiba analisá-la em regime estático. Esta 
asserção mostra bem a interdependência entre 
os dois tipos de análise estrutural e é facilmente 
verificada através da leitura dos trabalhos apre- 
sentados em reuniões internacionais sobre cál- 
culo automático de estruturas [22 a 28). 

Como exemplo de aplicações recentes do mé- 
todo de análise modal refere-se o cálculo sis- 
mico de barragens abóbada, em que a determi- 
nação da matriz deformabilidade foi feita por 
via analítica [29, 30] e experimental [31], de 
barragens de terra [32] e de estruturas de edifí- 
cios de painéis préfabricados [33]. 

Indica-se na fig. 8 os valores das frequências 
próprias e modos de vibração de um edifício de 
6 pisos daquele tipo, considerando-o encastrado 
e assente sobre um meio elástico. 

No que diz respeito à análise dinâmica de 
estruturas em regime não linear, a aplicação dos 
métodos de integração «passo a passo» para O 
estudo do comportamento de osciladores de um 


a) Edifício encastrado 


T,:0,0835 
f4 = 124 Hz 


52 Hz 


08 Hz 


; 


(ed 


e vários graus de liberdade conduz a um volume 
de cálculo muito superior ao do estudo do com- 
portamento em regime linear, pelo que se torna 
imprescindível o emprego de computadores, 

De entre os vários campos de análise diná- 
mica de estruturas em que o cálculo automático 
tem desempenhado um papel importante, salien- 
ta-se o da Engenharia Sísmica, cujo enorme 
progresso na última década se deve em grande 
parte ao desenvolvimento dos computadores. 
Basta folhear as actas das Conferências Mundiais 
de Engenharia Sismica para verificar que a rea- 
lização da maioria dos trabalhos apresentados 
não teria sido possível sem o concurso de com- 
putadores, digitais ou analógicos [34, 35, 36). 

Apresenta-se na fig. 9 o diagrama desloca- 
mentos-tempo de dois osciladores bi-lineares de 
um grau de liberdade com amortecimento vis- 
coso, sujeito à acção de um sismo de forte inten- 
sidade (componente NS do sismo de El Centro, 
Califórnia, de 18 de Maio de 1940). A integra- 
ção passo-a-passo foi realizada no computador 
Elliott 803-B com que está equipado o Centro 
de Cálculo do LNEC. 


154 Hz 


191 Hz 212 Hz 


1º modo 6º modo 
b) Edificio apoiado sobre um meio elástico 
7,=:0,357s 
f -2.8Hz 12.6 Hz 78Hz 138 Hz 184 Hz 214 Hz 


vw bh 4a 


ie modo 


ONCE: 


6º modo 


Fig. 8 — Frequências próprias e modos de vibração de um edifício de painéis préfabricados de 6 pisos [33] 
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ACELEROGRAMA DO SISMO DE EL CENTRO, 18-5-1940 (NS) 


300 


3 


100 


ACELERAÇÃO, á(t), cm/seg? 
S 


DESLOCAMENTO RELATIVO x/xo 
1 


— Oscilador bi-linear: wo= 2H7;1:0,05;Xo:2.5em;r =0,12 
8 [= = Oscilador elasto-plástico: Wo=2H2;1:0,05;X0=2,5cm;r=O 


OSCILADOR DE UM 
GRAU DE LIBERDADE 
x(t) 


Wo * Frequência própria 
1 = Amortecimento relativo 


ra +98 
tg 


so TSE 


TEMPO, seg. 


Fig. 9 — Resposta de osciladores bi-lineares sujeitos à acção do sismo de El Centro de 1940 (NS) 


6 — CONCLUSÕES 


A análise de estruturas actuadas por solicita- 
ções em regime dinâmico sofreu nos últimos 
vinte anos um desenvolvimento considerável 
devido ao progresso dos conhecimentos nos 
vários domínios a que está ligada. 

A formulação dos vários métodos de que dis- 
põe é em geral muito simplificada quando se 
adopta a representação matricial. 

O desenvolvimento dos computadores electró- 
nicos exerce uma influência dupla sobre a aná- 
lise dinâmica de estruturas; por um lado torna 
viável a aplicação numérica dos métodos gerais 
da Teoria das Vibrações, permitindo resolver 
problemas de dimensionamento ; por outro lado 
impulsiona a investigação dos vários campos em 
que a análise dinâmica se pratica, particular- 
mente dentro da Engenharia Sísmica. 
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DISTÂNCIA MÍNIMA NOS CÓDIGOS DE BOSE-CHAUDHURI- 
-HOCQUENGHEM. MÉTODO PARA A DETERMINAÇÃO 
DE UM LIMIAR 


RESUMO 


Expõe-se a teoria que determina um método expedito 
para a determinação de um limiar para a distância mi- 
nima nos códigos de Bose-Chaudhuri. Utiliza-se para 
tanto um polinômio g, (X) cujo grau determina direc- 
tamente o limiar da distância minima. Apresenta-se a 
terminar um exemplo para se julgar da utilidade do 
método. 


Introdução 


Seja n um número natural ímpar superior à 
unidade e seja V um espaço vectorial linear com 
elementos 

aÃ, Vo» ma Cao) (1) 
que são vectores baseados no corpo de Galois 
GF (2). Define-se adição de vectores em V, 
sendo o vector soma o vector cujas componentes 
são a soma, em GF(2), das componentes dos 
vectores parcelas. 

Pode-se considerar qualquer código binário, A 
em que a dimensão das palavras é igual an, 
como um sub-espaço de V: 


A CV 


Consideremos a, uma palavra qualquer do 
código A ese a =(a,,a,,...,a,.,) impli- 
car que a, =(a, ,,a,,...,a, 3) também é 
uma palavra do mesmo código A, então este 
código é um código cíclico. 

É nos códigos de Bose-Chaudhuri, que são 
códigos cíclicos que se aplica a teoria exposta 
neste artigo. 

Vamos então considerar uma série de resulta- 
dos, com os quais seja possível um método 
expedito para a determinação do limiar da dis- 
tância mínima ou distância de Hamming nos 
códigos de Bose-Chaudhuri. 
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SYNOPSIS 


In this paper we present the theory used in a method 
for the determination of minimum distance bounds for 
Bose-Chaudhuri codes We define a polynomial g, (X) 
the degree of which is directly related to a minimum 
distance bound. At the end the presentation of an exam- 
ple shows the possibilíties of the method. 


Lema 1: 


Consideremos uma sequência em GF(2) a, a,, 

«y à, -y Obtida a partir dos valores iniciais 

ay à, ... ay, e utilizando a fórmula de recor- 
rência (ou equação às diferenças) 


acutb agm+...+ ba=o0 (2) 


em quei=-0,1,...,(n-k-1) e os coeficientes b, 
são elementos de GF (2). Se o polinómio 


E)=X Hb XHIF...+b (3) 


não tiver raízes múltiplas e dividir (1 + X") en- 
tão existem elementos univocamente determi- 
nados num corpo de extensão K tais que 


a=c W'+toWi+...ta Wi (4) 
em que W, W,..., W, são as raízes de f(X). 


Demonstração : 


k 
Sea,= 2 c Wo', então substituindo na fór 
mula de recorrência (2) vem 
k k 
2 q WtH+ b 2 Wçtii+... 
6=—1 


(1=1 


k 
+b, 2 qWi=o 
0=1 


| tar 


O wo (Wo + by Wjt-! + 
+. b)=0 (5) 


fi 
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Logo é necessário que 
Wa +b, Was +. . Hb =0 


para qualquer ?=1,2,...,k,istoé, Wi, W,..,W, 
são raízes do polinómio f(X) =Xk + b, Xk=! + 
+... b, e sendo em número de k, são todas 
as raizes. 

Como quaisquer k equações se podem utilizar 
para a determinação dos c,'s, vamos usar 


a=cWitoWi+...+taW: (6) 


com i = 0,1,2,..., (k—1). Na solução do sis- 
tema assim obtido recorre-se ao seguinte deter- 
minante 


W MW; W; e W 
| l 1 l 
WoW; W. e W, 
(7) 
WoW Wi; W. 
k—l k—l k—l k—l 
W) Wi w' Wi 


O determinante acima é um determinante de 
Van der Monde. Se os valores W,, Wo, Wy,..., 
W., que são as raízes de f(X), forem distintos 
o determinante de Van der Monde não é nulo e 
portanto os c;'s obtêm-se por uma transformação 
linear não-singular a partir de quaisquer k a;'s » 
isto é, univocamente. QED., 


Lema 2: 


O sub-espaço A do espaço vectorial V defi- 
nido através da fórmula de recorrência 


ces ja:a=(aç,a,,..., 20 P 


ax TD, ar + ces bd =0 
parai=0,1,... Eine JE 


é um código cíclico, cujo polinómio gerador é 


80X) = (1+X)/£ 09. 


Demonstração : 


Demonstremos em primeiro lugar que o código 
é cíclico. 
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Os W são todos raizes de índice n da unidade 
pelo facto de 14X” ser múltiplo de f(X) e 
portanto Wº = 1. 

Logo utilizando a expressão (3) 


k k 
, n i q i id 
As 2 aW”! = 2 MW, = 
0=1 0=1 
* . 
= É em 
= é CG W= a (8) 
,==1 
Verificámos que a, .; =a, isto é, os símbolos 


repetem-se ciclicamente com período n. 

Vamos agora mostrar que o polinômio gerador 
é g (X). Comecemos por ver que a fórmula de 
recorrência (2) 


| a +tbag,t...+tb,a=o0 
astrba, +t...+tb, a =0 
| dica E Da E cs + É Dra = 
corresponde ao produto do vector a ==(a,, a,, 


do,» ..«, à, 1) pela matriz de m linhas por 
(n — k) colunas 


fi b, 0 0 say ças A “i 
bia Bb dO 0 
Dig by b, 0 
o db MB o (9) 
1 b b, 0 
O 1 b, O 
0 0 0 se use 
O) O (O) a cos A 


Portanto o vector à é ortogonal a cada uma 
das colunas da matriz (9) e a conclusão é de 
que as palavras do código são ortogonais a 
(n —k) vectores que formam a base vectorial do 
código ortogonal A*. Qualquer vector a* e A* 
é ortogonal aos vectores a € A e vice-versa, 

Utilizando polinómios módulo 1X” para 
simbolizar vectores, a base do código ortogo- 
nal A* é formada pelos seguintes polinómios 
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* th x Eb xt t...tbp= 


= f (X) 
XE de by xe +b, ipa do co dp bpX= 
= X f(X) 
Et 2 gy SP vero HÃO 
= X2 £ (X) 
XxX +bhX2 +... + XT = 
X n-k—l E (X) 


Como o polinómio f(X) é de grau k os n-k 
polinómios formam a base do código ortogonal 
A* e f(X), sendo o polinómio gerador de menor 
grau da base, é denominado por polinómio gera- 
dor de A“, Finalmente como os polinómios gera- 
dores dos códigos cíclicos ortogonais A e A* de 
dimensão n são caracterizados pelo facto do seu 
produto ser igual a 1 + X” e como o polinómio 
gerador de A* é f (X), então o polinómio gera- 
dor de A será (1 + Xº)/ f(X), que é exacta- 
mente o polinómio g(X). QED. 

A definição do código cíclico utilizando a fór- 
mula de recorrência indica-nos que os primeiros 
k símbolos são os símbolos portadores de infor- 
mação e determinam completamente os restantes 
(n-k) símbolos que constituem a parte redun- 
dante das palavras. 


Definição 1: 


Se Z; é uma raiz primitiva de índice n da uni- 


dade , o polinómio f(X) fica completamente 
determinado por um conjunto de números in- 
teiros 


E (fZ) = |e:o <e<n,f(Z*) =0 | = 


= | Cp Cos o os Ce! 


em que o número de elementos é igual ao grau 
de f(X). 

Definição 2: 

Vamos escolher uma raiz primitiva de índice 
n da unidade Z tal que emax, O maior inteiro em 


E (f,Z), seja quando muito igual ao maior inteiro 
em E (f,Z,), sendo Z, qualquer outra raiz primi- 


tiva de índice n da unidade. 
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É utilizando estas definições que vamos con- 
tinuar. 


Lema 3: 


Dada uma raiz primitiva de índice n da uni- 


dade, Z,, para cada palavra do código A, a = 
= (a, ,à,,...,ã, 1), existe um único polinó- 
mio g, (X) com coeficientes num corpo de exten- 


são K tal que 
a=s, (Zi) para j=0,1,2,..., (n—1). 
O grau do polinómio g, (X) é menor ou igual 
ao maior inteiro em E (f, Z,). 
Demonstração : 


Consideremos a raiz de índice n da unidade Z,. 
Então as raízes do polimónio f (X) são 


ZÃ, Msc E 


e não há nenhuma repetida visto f (X) não ter 
raízes múltiplas, 
Segundo o Lema 1 existem coeficientes c,, C,, 


« y Cy tais que 


k k 
a= 2 a(Z)i= 2 a(Z)º (1) 
1=1 6=1 
e portanto o polinómio 
k 
By VÃ) = 2 qX* (12) 
0= 1 


é univocamente determinado em face de Z,, uma 
vez que os coeficientes c,, €,, . . . , C, São únicos 
para qualquer palavra a e A (Lema 1). 
Deve-se notar que g, (X) éiguala 10ouao 
quando X é substituído por qualquer das raízes 


primitivas de índice n da unidade. 
k 
Da expressão g, (X) = 2 
1=1 
que o grau de g, (X) é menor ou igual ao maior 
dos inteiros em E (f, Z,)). QED. 
Se a raiz primitiva de índice n da unidade 
utilizada for Z então o grau de g, (X) é menor 


ou igual a €qav 


e . 
c, X% se conclui 
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Definição 3: 


Se g(X)=(1+Xº)/f(X) tiver entre as suas 
raízes YU, Yº, Yº,..., ySom 
tiver Yºº então o código cíclico obtido através 
da fórmula de recorrência baseada em f(X), é 
um código de Bose-Chaudhuri associado a Y e 
ao número d,. 


6 ia 
mas já não con- 


Apresentemos agora o resultado mais impor- 
tante sob a forma de teorema. 


Teorema 1: 


A distância de Hamming para o código de 
Bose-Chaudhuri associado a Z”! e ao número d, é 


w(A)2d,=n— €a 


se Z for a raiz primitiva de índice n da unidade 
indicada na Definição 2. 


Demonstração : 


Considerando 
w (a) = [ número de símbolos 1 na palavra a € A ) 


então a distância de Hamming é o menor valor 
que w (a) toma em A, com exclusão de a = 0 = 
= (0,0,...,0). Se representarmos a distância 
de Hamming para o código A por w(A), então 


w(A) =min(número de símbolos 1 na palavra a! = 
acÃ 


a+0 


= min (n) fnúmero de símbolos O na palavra a!) = 
aca 
a+0 


= n — max (número de símbolos O na palavra 8! — 
ach 
a+ 0 


= n — max (número de raízes de g, (X) em Z) 
acÃ 
a0 


sendo Z o grupo multiplicativo das raizes de 
índice n da unidade. 

Utilizando para definir g, (X) uma raiz qual- 
quer de índice n da unidade, Z, o número de 
raízes de g, (X) será no máximo o maior inteiro 
no conjunto E (f,Z,); seja m, o maior inteiro em 
E (f, Z,). Interessa, por conseguinte, o menor dos 


4 
m, que é exactamente o valor e, ..- 


Logo 
max f número de raízes de g, (X) em Z | < emas 
ac À 
a+0 


Como (Z-!)dº não é raiz de g (X), terá de ser 
raiz de f(X). Então d, é o inteiro de menor va- 


lor absoluto em E (EZ O. Mas como 
(27! to 74 = 7"-4% 


é uma raiz do polinómio f (X), então (n—ds) não 
só pertence a E (f,Z) como é o inteiro de maior 
valor absoluto; daí a conclusão 
max TT A — d, 
ou, finalmente, 
w (A) > d, = N — max 
QED. 


Corolário : 


Para o código de Bose-Chaudhuri associado a 
Y e ao número d tem-se 
ds n— max 


Demonstração 


Sendo W = de então E (£,W) contém como 
inteiro de maior valor absoluto (n—d) e sendo 


A ce n—d vem 


dA n= Cm 


QED. 


Exemplo de aplicação destes resultados na 
determinação do limiar da distância de Hamming 
para alguns códigos de Bose-Chaudhuri com 
n = 21. 

21 : 

Factorizando 1 + Xº pelo método dos poli- 


4 . . , . (2) 
nómios ciclotómicos 


1 + XE = G, (X). G; (X). G; (2X). G; (X). 6, 09). 
AB; 00) 
em que o índice indica o grau do polinomio. 
Então 
E(G, Y)=[1,2, 4,8, 16,11) 


E(G., Y)=[5, 10, 20, 19, 17, 13) 
E(G, YW)=(3, 6,12] 
E(G;,Y)=(9, 18, 15] 

| 
E(G, YW=(7, 14) 
E(G: W=(0] 
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Fazendo agora, por exemplo, 

a) FX) =G(X). G;(X). G,(X). G, (X) 
obtém-se e, = 16 e portanto d,=5; 

b) F(X)=G. (X). G,(X). G, (X) 
obtém-se e ay = 14 e portanto d, = 7; 

o) E(X) = G,(X). G, (X) 


obtém-se eax = 12 e portanto d, = 9. 


Os códigos correspondentes a estas três alíneas 
são, respectivamente, (21, 12), (21, 6) e (21, 4), 
Muitas outras combinações de factores se pode- 
riam tentar para formar f(X), mas os códigos 
obtidos não seriam tão eficientes e portanto não 
têm qualquer interesse prático. 


Conclusões 


Este artigo é o ínicio de um trabalho sobre 
códigos de Bose-Chaudhuri dedicado à determi- 
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nação e estabelecimento de limiares na distância 
de Hamming. Na sequência deste artigo e utili- 
zando o material introduzido por ele, tem-se 
procurado alguns casos particulares em que o 
limiar de w(A) se pode estabelecer com mais 
rigor, isto é, ir além de afirmar simplesmente 


que excede (n — e...) - Tal foi o assunto da comu- 


nicação apresentada no Simpósio Internacional 
de Teoria da Informação, realizado em Itália 
(San Remo), em Setembro de 1967, com o título 
«On a Class of Wide-Sense Binary BCH Codes 
Whose Minimum Distances Exceed the BCH 
Bound». Prevê-se para um futuro próximo a 
publicação, em português, do texto desta comu- 
nicação. 


BIBLIOGRAFIA 


(1) A new treatment of Bose-Chaudhuri Codes — H.F. 
Mattson, G. Solomon — Journal of Society for Indus- 
trial and Applied Mathematics, SIAM, vol. 9, n. 4, 
p. 654-669 (Dez. 1961); 

2) Error-Correcting Codes — W.W. Peterson — MIT 
Press, Cambridge, Massachusetts, 1961. 


217 


V SEMANA DE ENGENHARIA CIVII 
9-18 DE MARÇO DE 1968 


PROGRAMA 
9 de Março 


11 horas — Sessão de Abertura no Salão Nobre do 1.5.T. 


|| de Março 


11 horas — Salão Nobre do 1.5.T. — Conferência sobre «A Noção de Vão Limite e a Sua Uti- 
lização no Projecto de Estruturas», pelo Prof. Armando Campos e Matos, da 
Faculdade de Engenharia do Porto. 


|2 de Março 


11 horas — Salão Nobre do 1.5.T. — Conferência sobre «Segurança nos Estaleiros», pelo Eng.º 
Pinto de Almeida, da H.I.C.A. 

16 horas — Conferência sobre «Time-Sharing», pelo Dr. Barreto Braga, consultor científico 
da Bull. 


13 de Março 
10 horas — Visita ao Centro Mecanográfico do Ministério das Finanças. 


14 de Março 
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Equipamentos electromecânicos completos para centrais e fábricas 
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